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ABSTRACT 


Presented  in  this  paper  is  a  new  matrix  representation  of  classical  electromagnetic 
theory.  The  basis  of  this  representation  is  a  space-time,  eight-by-eight  differential  matrix 
operator.  This  matrix  operator  is  initially  formulated  from  the  differential  form  of  the 
Maxwell  field  equations  in  vacuum.  The  resulting  matrix  formulation  of  Maxwell's 
equations  allows  simple  and  direct  derivation  of  the  electromagnetic  wave  and  charge 
continuity  equations,  the  Lorentz  conditions  and  definition  of  the  electromagnetic 
potentials,  the  Lorentz  and  Coulomb  gauges,  the  electromagnetic  potential  wave 
equations,  and  Poynting's  conservation  of  energy  theorem.  A  four-dimensional  Fourier 
transform  of  the  matrix  equations  casts  them  into  an  eight-dimensional  transfer  theorem. 
The  transfer  function  has  an  inverse,  and  this  allows  the  equations  to  be  inverted.  This 
inversion  expresses  the  fields  directly  in  terms  of  the  charge  and  current  source 
distributions,  i.e.,  without  the  need  for  calculating  intermediary  potentials.  This 
inversion  formula  is  new,  for  the  general  scenario  of  space-  and  time-dependent  sources. 
A  simple  pedagogical  example  is  included  illustrating  use  of  the  formulation. 
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1.0  INTRODUCTION 


Present  day  mathematical  descriptions  and  computations  in  the  field  of  classical 
electrodynamics  rely  heavily  on  the  use  of  vector  calculus.  Matrix  calculus,  on  the  other 
hand,  finds  widespread  use  in  other  areas  of  physics,  engineering,  and  applied 
mathematics.  Structural  analysis,  vibrational  analysis,  electrical  circuit  analysis,  and  the 
theory  of  elasticity  are  just  a  few  of  the  disciplines  (Pipes,  1963)  where  matrix  calculus 
has  been  successfully  employed  routinely.  This  paper  presents  a  new  formulation  of 
electromagnetic  theory  from  a  matrix  calculus  point  of  view.  We  will  show  the 
formulation  arises  from  use  of  a  skew-Hermitian  space-time  eight-by-eight  differential 
matrix  operator 
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By  "skew-Hermitian"  we  mean  [M]  H  =  -[M]  ,  where  [M]  H  denotes  the  Hermitian 
(complex  conjugate)  transpose  of  [M].  The  notation 

x  =  id,  i  =  I  (2) 


is  adopted  for  describing  the  temporal  component  of  a  space-time  point.  All  basic 
physical  quantities  cf  the  theory  will  be  four-dimensional,  with  a  complex  fourth 
component  as  in  equation  2. 
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2.0  DERIVATION  OF  THE  [M]  MATRIX 
REPRESENTATION 


The  fundamental  equations  of  classical  electromagnetic  phenomena,  namely  the 
Maxwell  field  equations,  serve  as  our  starting  point.  In  this  section  the  space-time 
operator  [M],  which  is  of  paramount  importance  to  the  remainder  of  the  paper,  is 
constructed  from  the  Maxwell  field  equations.  In  the  Gaussian  system  of  units,  the  four 
Maxwell  field  equations  in  vector  form  are  given  by  the  following  (Jackson,  1962): 


Ampere-Maxwell  law 

VxH(r.f)  =  (-)|jD(r,0  +  (— )f(r,t)  (3a) 

c  dr  c 

Gauss'  law  for  electricity 

V*D  (r,  t)  =  (4k)  pe  (r,  t)  (3b) 

Faraday's  law 

VxE(r.l)  =-  (i)  §-B  (r,  I)  -  (-)f' (r,  t)  (3c) 

C  of  C 

Gauss'  law  for  magnetism 

V«B(r,t)  =  (4K)pm(r,f)  (3d) 


The  order  in  which  the  four  equations  appear  above  is  important  in  the  construction  of 
the  space-time  operator  [M],  In  rectangular  coordinates 
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JL  A) 
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(4) 


and 

r  =  ( x,y,z ) . 


(5) 
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The  physical  quantities  appearing  in  the  Maxwell  field  equations  are 


£  (r,  t)  -  Electric  field 

D(r,t )  -  Electric  displacement 

B  ( r ,  t )  -  Magnetic  induction 

H  (r,  t)  -  Magnetic  field 

f  ( r,t )  -  Electric  current  density 

Jm  (r,t)  -  Magnetic  current  density 

pe  (r,  t)  -  Electric  charge  density 

pm  ( r ,  t )  -  Magnetic  charge  density 

V  -  Del  operator 

V«  -  Divergence  operator 

Vx  -  Curl  operator 

(r,  t)  ~  Space-time  point 

c  -  Speed  of  light  in  vacuum 


Both  magnetic  charge  and  current  densities  (Magid,  1972)  have  been  included  in 
Maxwell's  equations  for  purposes  of  completeness.  These  densities,  of  course,  may  be  set 
equal  to  zero  since  magnetic  monopoles  have  not  yet  been  discovered  in  nature. 

The  electric  displacement  and  electric  field,  as  well  as  the  magnetic  induction  and 
magnetic  field,  are  related  (Jackson,  1962)  through  the  expressions 

D  (r,  t)  =  E  (r,  t)  +  (4n)  P  (r,  t)  (6a) 


B(r,t)  =  H(r,t)  +  (4n)M(r,t)  (6b) 

where 

P(r,t)  -  Macroscopic  polarization 
M(r,t)  -  Macroscopic  magnetization. 


In  this  paper,  dielectric  and  magnetic  materials  are  not  considered.  Hence,  we  let 


P  (r,  t)  =  0  (7a) 

M  (r,  t )  =  0.  (7b) 


A  subsequent  paper  will  treat  dielectric  and  magnetic  materials  where  both  the 
macroscopic  polarization  and  magnetization  vectors  are  non-zero. 

With  the  use  of  equations  6  and  7,  the  four  Maxwell  equations  for  vacuum  can  be 
rewritten  as  eight  scalar  equations.  That  is. 


Ampere-Maxwell  law 
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Gauss'  law  for  electricity 

J$-E  +J$-E  +|-E  =  (  —  )cpe 
dx  x  dy  y  dz  z  v  c  y 

Faraday's  law 

5y£2"K£y  y  aiH*  - 
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Gauss'  law  for  magnetism 

|-H  +|~H  +|-H  =  (  —  )cpm. 

dx  x  dy  y  dz  z  c 


(8a) 
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The  eight  scalar  equations  can  now  be  combined  into  a  single  matrix  equation 
involving  an  eight-by-six  differential  matrix  operator: 
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Other  authors  have  formed  electromagnetic  matrix  operators  similar  to  equation  9,  but 
either  under  less  general  conditions,  or  with  different  goals  in  mind  (Mayes,  1967;  Eisele 
and  Mason,  1970).  Mayes  (1967)  assumes  the  special  case  of  time-harmonic  fields  due  to 
radiating  antennas,  so  that  time  is  suppressed  in  the  formulation.  The  effect  is 
suppression  of  Maxwell's  equations  3b  and  3d,  resulting  in  a  six-by-six  matrix  operator 
formulation  instead  of  eight-by-six  matrix  operator  formulation  (equation  9).  Eisele  and 
Mason  (1970)  assume  zero  source  strengths,  and  form  complex  field  quantities, 
Ex  +  iE  ,  Ez,  Hx  +  iHy,  Hz  ,  with  the  goal  of  showing  Maxwell's  equations  reduce  to  a 
Dirac-like  equation  for  a  four-component  electromagnetic  spinor  field  with  zero  mass. 

Also,  no  other  authors,  as  far  as  we  know,  take  the  step  described  next  of  broadening 
the  long  (eight-by-six)  matrix  in  equation  9  into  an  associated  square  (eight-by-eight) 
matrix,  with  attendant  theoretical  benefits. 

We  are  now  in  a  position  to  construct  the  desired  eight-by-eight  space-time  operator 
[M].  The  key  to  this  construction  is  to  first  replace  the  eight-by-six  matrix  operator 
appearing  in  equation  9  by  an  eight-by-eight  matrix  operator  containing  two  additional 
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columns.  That  is. 


'  a 

3x 

0 

0 

0 

a 

az 

A 

™18 

lEv 

f 

Jx 

0 

a 

ax 

0 

™24 

a 

dz 

0 

a 

3x 

™28 

X 

iE 

f 

0 

0 

a 

ax 

™34 

_d_ 

_a_ 

a^ 

0 

™38 

y 

iE 

z 

y 

fz 

a_ 

dx 

a_ 

A 

dz 

m44 

0 

0 

0 

™48 

0 

=  <y> 

icpe 

0 

_d_ 

az 

A 

™54 

_a_ 

dx 

0 

0 

™58 

H 

X 

-if 

Jx 

a 

az 

0 

a 

dx 

™64 

0 

a 

ax 

0 

™68 

H 

y 

A 

_a_ 

dy 

A 

dx 

0 

rh14 

0 

0 

_a_ 

ax 

m78 

H 

z 

< 

0 

0 

0 

m84 

A 

3x 

A 

a_ 

3z 

™88 

spm_ 

(10) 


Columns  4  and  8  are  the  newly  inserted  columns.  Their  elements  rh-j  are  as  yet 
undefined.  To  null-out  the  effects  of  the  two  new  columns,  additional  zeros  must  be 
inserted  in  the  electric  and  magnetic  field  vector  at  positions  4  and  8.  This  eight-by-eight 
matrix  operator  representation  of  Maxwell's  equations  is  now  entirely  equivalent  to  the 
original  vector  representation,  regardless  of  choice  of  the  added  elements  m-j  .  We  next 
choose  useful,  specific  values  for  the  added  m-j  . 

The  motivation  for  choice  of  the  added  elements  m-  is  as  follows.  The  off-diagonal 
known  elements  in  the  eight-by-eight  matrix  have  an  antisymmetric  arrangement.  For 
example,  element  m  5  is  and  m52  is  .  This  arrangement  suggests  that  we 
choose  the  added  elements  m  ■  ■  such  that  antisymmetry  is  preserved.  This  preservation 
requires  that  we  choose  the  off  ^diagonal  elements  m  -  as  follows: 
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and 
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(Actually,  assigning  =  m48  =  0  is  arbitrary;  these  elements  can  be  assigned  as 
m84  =  b  and  m48  =  -b  ,  with  b  arbitrary.  We  choose  b  =  0  for  simplicity.) 

Next,  the  diagonal  elements  rh ^  and  m88  have  to  be  specified.  If  these  elements 
are  chosen  according  to  the  equation 

mu  =  m88  =  -J^  (13) 

then  the  conditions  for  the  Lorentz  gauge  are  satisfied.  This  will  be  seen  in  the  next 
section.  If,  however,  these  elements  are  chosen  according  to  the  equation 

"*44  =  "*88  =  0  (14) 

then  the  conditions  of  the  Coulomb  gauge  are  satisfied.  See  appendix  A  for  more  details 
on  the  Coulomb  gauge.  Substituting  the  Lorentz-gauge  choice,  equation  13,  and  the 
antisymmetric  choices,  equations  11  and  12,  in  the  general  form,  equation  10,  leads  to  the 
following  eight-by-eight  matrix  representation  of  Maxwell's  equations: 
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The  left-hand  matrix  is  the  matrix  in  equation  1  we  previously  called  [M].  Derivation 
of  this  matrix  form  (equation  15)  of  the  Maxwell  field  equations  is  the  chief  result  of  this 
section.  This  matrix  will  be  used  throughout  the  remainder  of  the  paper.  This  matrix 
permits  simple  matrix  operations  upon  equation  15  to  derive  the  major  phenomena  of 
electromagnetic  theory.  This  derivation  is  the  subject  of  section  3.  It  also  permits  direct 
inversion  for  the  fields  £  and  H  in  terms  of  the  sources,  without  the  need  for 
calculating  intermediary  potentials.  This  inversion  is  the  subject  of  sections  4  and  5. 

Equation  15  can  be  rewritten  in  a  more  compact  form 
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In  summary,  the  matrix  equation  16  is  equivalent  to  the  Maxwell  field  equations  in 
vacuum.  Appendix  B  shows  the  matrix  representation  (equation  16)  is  in  covariant  form. 
This  development  will  continue  in  section  4  where  four-dimensional  Fourier  transforms 
will  be  taken  in  order  to  achieve  inversion  of  the  matrix  field  equations. 

3.0  MATRIX  FORMULATION  OF  OTHER 
ELECTROMAGNETIC  PHENOMENA 

The  matrix  representation  (equation  16)  of  the  Maxwell  field  equations  will  imply,  in 
simple  fashion,  some  milestone  effects  of  electromagnetic  theory.  Simple  matrix 
operations  upon  the  space-time  operator  [M]  constructed  will  be  shown  to  derive  (a)  the 
electromagnetic  wave  and  charge  continuity  equations,  (b)  the  Lorentz  conditions  and 
definition  of  the  electromagnetic  potentials,  (c)  the  electromagnetic  potential  wave 
equations,  and  (d)  Poynting's  conservation  of  energy  theorem. 

3.1  ELECTROMAGNETIC  FIELD  WAVE  EQUATIONS  AND 
CHARGE  CONTINUITY  EQUATIONS 

Multiply  both  sides  of  matrix  equation  16  by  the  complex  conjugate  of  the  space-  time 
operator  [M].  This  gives 
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and  □ 2  is  the  D'Alembertian  operator  (Feynman,  Leighton,  and  Sands,  1964)  defined  by 
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The  matrix  equation  20  is  equivalent  to  eight  scalar  equations.  These  scalar  equations  are 
equivalent  to  the  vector  electromagnetic  field  wave  equations 
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and  the  electric  and  magnetic  charge  continuity  equations 

V.f  +  jLpe  =  0  (24a) 

of 

v#/m+itpm  =  °*  (24b) 

In  summary,  a  single  matrix  multiplication  of  both  sides  of  the  Maxwell  field 
equations  (equation  16)  by  the  complex  conjugate  of  the  space-time  operator  [M]  yields  a 
single  matrix  equation  which  is  equivalent  to  both  the  electromagnetic  field  wave 
equations  and  the  charge  continuity  equations. 

3.2  ELECTROMAGNETIC  POTENTIALS  AND  LORENTZ 
CONDITIONS 


We  found  before  that  the  space-time  operator  [M]  defines  the  eight  scalar  Maxwell 
field  equations.  We  now  observe  that  the  complex  conjugate  of  [M]  also  provides  the 
definition  of  the  electromagnetic  fields  in  terms  of  the  familiar  vector  and  scalar 
potentials.  In  particular,  for  the  Lorentz  gauge,  the  matrix  equation 
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is  equivalent  to  eight  scalar  equations.  Six  of  these  are  equivalent  to  the  two  vector 
equations  which  define  the  relationship  between  the  electromagnetic  potentials  and  the 


electromagnetic  fields 

E  _  _V(pe- (-)$-Ae-VxAm  (27a) 

c  at 

H  =  -V<pm-  (-)~-Am  +  VxAe.  (27b) 

c  at 

The  remaining  two  scalar  equations  correspond  to  the  Lorentz  conditions 

V.Ae+  (i)^e  =  0  (28a) 

V.Am+  (-)^m  =  0.  (28b) 

c  at 


Results,  equations  27  and  28,  are  easily  verified  by  explicit  multiplication  of  the  right- 
hand  side  of  matrix  equation  25. 

It  is  interesting  that  the  two  vectors  [p j] ,  [p2]  are,  by  their  definitions  (equation 
26),  relativistic  four-vectors.  A  comparison  of  matrix  equations  16  and  25  discloses  the 
curious  fact  that  the  same  operator  [M]  that  operates  upon  the  fields  f =1,2  to  give  the 
sources  s .  also  operates  upon  the  potentials  p-  to  give  the  fields.  There  is  an 
interesting  symmetry  between  the  two  effects. 


3.3  ELECTROMAGNETIC  POTENTIAL  WAVE  EQUATIONS 


We  know  the  electromagnetic  vector  and  scalar  potentials  satisfy  inhomogeneous 
wave  equations.  Once  again,  we  show  this  by  simple  matrix  multiplication.  First 
substitute  the  matrix  expression  (equation  25)  into  the  matrix  representation  (equation 
16)  of  Maxwell's  equations.  This  gives 


* 


M,  M2 

M1  M2 

Pi 

M2  M1 

m2  m1 

Vi 

(29) 


However,  for  the  Lorentz-gauge  choice  of  [Mj, 


1 

"1 

* 

_  _ 

* 

- 

M1  M2 

Mj  M2 

Mj  /Vi 2 

Mj  M, 

D  0 

M2Ml 

M2Mj 

M2  Mj 

M2  Mj 

0  D 

(30) 
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by  direct  evaluation  of  the  matrix  product.  Then,  by  equations  29  and  30,  the  following 
matrix  representation  of  the  electromagnetic  potential  wave  equations  is  obtained 


D  O' 

P\ 

0  D 

Pi 

(31) 


By  explicit  evaluation  of  the  left-hand  matrix  product,  this  single  matrix  equation  is 
equivalent  to  eight  scalar  equations.  Six  of  these  scalar  equations  are  equivalent  to  the 
following  pair  of  inhomogeneous  wave  equations  for  the  electric  and  magnetic  vector 
potentials 


V2Ae-(l)-^Ae  =  -(— )f 
c2  dt2  c 


(32a) 


V2Am-(4) 


,m 


dt' 


=  -A/m 

c 


(32b) 


The  remaining  two  represent  the  inhomogeneous  wave  equations  for  the  electric  and 
magnetic  scalar  potentials 


vV-(4>-^>e  =  -(4*)pC 

cL  dr 

VW"—  (i)-2V  =  -(4„)pm 

C  dt 


(33a) 

(33b) 


3.4  ENERGY  CONSERVATION;  POYNTING'S  THEOREM 

The  law  of  conservation  of  energy,  often  called  Poynting's  theorem  (Jackson,  1962),  is 
an  important  milestone  of  electromagnetic  theory.  Again  simple  matrix  manipulation  of 
equation  16  will  accomplish  the  derivation.  Multiply  both  sides  of  equation  16  by  the 
Hermitian  conjugate  of  the  electromagnetic  field  vector,  fl  .  Directly, 


>1 


h 


Mi  M2 
M2Mj 


(34) 
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This  is  a  single  scalar  equation  representing  Poynting's  theorem 


V*S  +  $-u  +  E*f  +  H*Jm  =  0  (35) 

ot 

if  we  define  quantity  u  to  obey 

u=  (-L)  (  E  •  E  +  H.fl  ),  (36) 

where  u  is  the  total  energy  density  of  the  electromagnetic  fields,  and  define  S  by 

S-<£>  (  E  x  H  ),  (37) 

where  S  is  the  Poynting  vector  representing  energy  flow.  Again,  this  is  easily  verified 
by  explicit  multiplication  out  of  the  matrix  products  in  equation  34. 


4.0  ELECTROMAGNETIC  FIELD  FOURIER 
TRANSFORM  REPRESENTATION 


This  section  gives  a  method  for  inverting  the  Maxwell  field  equations.  The  Maxwell 
field  equations  have  been  written  in  a  compact  matrix  form  involving  the  square  eight- 
by-eight  space-time  differential  operator  [Ml.  With  the  use  of  four-dimensional  Fourier 
transform  theory,  the  matrix  form  of  Maxwell's  equations  can  be  easily  cast  from  the 
differential  form  into  an  algebraic  form.  The  algebraic  form  is  an  eight-dimensional 
transfer  theorem  in  K  -  space  ,  linearly  connecting  the  Fourier  transform  of  the  fields  to 
those  of  the  sources.  The  connection  is  through  an  eight-dimensional  transfer  function, 
whose  inverse  is  known.  This,  in  turn,  permits  inversion  of  the  field  equations.  The  fields 
are  then  directly  known  in  terms  of  the  charge  and  current  source  distributions.  (By 
"direct"  we  mean  without  the  need  for  calculating  intermediary  potential  functions.) 
Alternative  computation  of  the  electromagnetic  potentials  from  the  charge  and  current 
source  distributions  is  presented  in  section  7. 

A  three  dimensional  Fourier  transform  approach  to  direct  inversion  was  championed 
by  Mayes  (1985).  This  inversion  was  for  the  special  case  of  time-harmonic  fields  due  to 
antennas.  Also,  inversion  was  made  by  recourse  to  Green's  functions  (i.e.,  use  of  vector 
calculus)  rather  than  by  direct  use  of  matrix  transfer  theorem  (equation  45). 
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We  first  define  the  four-dimensional  Fourier  transform  pair 


(38a) 


(38b) 


where 

R=  (x,y,z,x)  =  (r,x)  and  K  =  (kx,k  ,kz,kx)  =  (k ,kx)  (39a) 

dR  =  dx  dy  dz  dt  and  dK  =  dk  dk  dk  da  (39b) 

^  y  ^ 

x  =  ict  and  kx  =  ica/c.  (39c) 

Next,  multiply  both  sides  of  equation  38b  by  the  space-time  operator  [M] 

dK.  W 


M,  M2 
M2Ml 


Mj  M 
M-  M 


+  i(JC.R) 


/i 

/2 


Operator  [M]  was  absorbed  under  the  integral  sign  since  the  integration  is  with  respect 
to  K ,  while  [M]  contains  space-  and  time-derivatives,  independent  of  K  .  We  can  easily 
show  that 


MtM2 

M2Mj 


e+i(K*R ) 


K  i  *2 
*2*1 


+  i(K»R) 


(41) 
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where 


*  1*2 
k2k, 


0  °  -kx  °  -kz  ky  0 

0  0  -*y  *2  0  ~kx  0 

0  0  -k%-kz-ky  k%  0  0 

*  k  k  -k,  0  0  0  0 

x  y  z  x 

0  -k  k  0  -JL  0  0  -it 

2  y  x  x 

k  0  -k  0  0  -k  0  -k 

z  x  x  y 

-it  it  0  0  0  0  -it  -it 

y  x  x  z 

0  0  0  0  kx  ky  kz  -k%  \ 


Then  by  equation  40, 


(42) 


M2 

/l 

M2M1 

h 

(-^) 

4k2 


m: 


K2 

k2k, 


f  1 

/2 


^KtR)dK. 


(43) 


Next  we  denote  the  Fourier  transform  of 
transpose.  Equation  16  becomes 


[SlS2l  by 


,  where  T 


is  the 


M,  M2 

h 

M2Mj 

f_2 

(44) 


Setting  equations  43  and  44  equal  to  one  another  gives  the  following  Fourier  transform 
representation  of  Maxwell's  equations 


>z 1 

* 

^  1 

fl 

k2  kx 

h. 

(45) 


This  is  an  electromagnetic  transfer  theorem,  the  direct  counterpart  to  well-known 
transfer  theorems  in  optics  (Goodman,  1968)  circuit  theory,  etc.  The  left-hand  side  matrix 
[K]  is  an  eight-by-eight  transfer  function.  It  has  an  inverse,  as  discussed  next. 


Transfer  function  matrix  [K]  has  elements  that  are  not  derivative  operators.  The 
elements  of  [K]  are  purely  algebraic,  i.e.  pure  numbers  in  specific  cases.  Therefore,  [K] 
can  potentially  have  an  inverse  in  contrast  to  matrix  [M]  in  direct-space  representation 
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(equation  16),  which  consists  of  differential  operators  and  hence  cannot  have  an  inverse 
in  the  ordinary  sense.  The  very  reason  for  taking  a  Fourier  approach  was,  in  fact,  to 
achieve  a  representation  with  a  realizable  inverse.  The  inverse  of  [K]  is  found  to  exist, 
and  is  given  by  (see  appendix  C) 


.-2 


Kl  K2 

K2Ki 


* 


(-- 2) 
*2 


Kl  K2 
*2*1 


(46) 


where 


*2  =  k*  +  k2+k2-  (a)2/c2). 


(47) 


This  is  easily  verified  by  explicitly  evaluating  the  matrix  product  of  the  right-hand  side 
of  equation  46  with  [K].  The  result  is  a  diagonal  identity  matrix,  as  required. 


With  [if]  1  known,  equation  45  may  be  inverted  for  the  fields  in  Fourier-space 


— 

-1 

♦ 

_ 

fl 

= 

A 

c 

*1  *2‘ 

Si 

-f-iVl 

*1  *2 

Si 

?2 

*2  *L 

h 

V  cK2) 

*2  *1 

S2 

by  identity  (equation  46).  Finally,  substituting  equation  48  into  equation  38b  gives  the 
fields  per  se, 


f\ 

h 


*1*2 

*2*1 


s2 


e+i«-R)dK. 


(49) 


The  poles  at  *  =  0  are  avoided  by  recourse  to  Cauchy  principal  value  integration 
(Jackson,  1962)  in  analytic  problems;  or  by  simply  avoiding  the  points  *  =  0  during 
numerical  integration  of  equation  49  in  nonanalytic  problems. 

Briefly  consider  experimental  cases  where  the  sources  [Sj]  ,[s2]  are  known  with 
inevitable  error,  e.g.,  due  to  random  noise.  In  this  case  equation  45  becomes  a  typical  ill- 
posed,  inverse  source  problem  (Ferwerda,  1978).  Inversion  formula  (equation  49)  would 
then  have  to  be  modified  so  as  to  avoid  unduly  magnifying  the  noise  at  points  *  ->  0. 
Appropriate  window  functions  that  smoothly  go  to  zero  as  *  — >  0  could  be  used,  e.g.,  to 
modulate  the  integrand  of  equation  49  for  this  purpose.  Or,  if  statistical  knowledge  of 
power  spectra  (\s-  (K)  |  ),  i=  1,2  is  forth  coming,  this  power  spectra  could  be  used  to 
form  a  Wiener  filter-based  inverse  solution. 
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We  have  thereby  achieved  our  goal  of  solving  for  the  electromagnetic  fields  directly 
in  terms  of  the  sources.  In  fact,  equation  49  shows  that  the  four-dimensional  Fourier 
transforms  of  the  sources  must  be  known  as  inputs.  In  cases  where  these  inputs  cannot  be 
computed  analytically,  numerical  use  of  the  four-dimensional  Fast  Fourier  transform 
(4D-FFT)  can  be  made.  Interestingly,  rectangular  coordinates  can  then  always  be  used, 
regardless  of  the  geometrical  shapes  of  the  sources. 

The  Fourier-based  approach  (equation  49)  to  inversion  appears  to  be  a  new  one  in  the 
general  scenario  of  space-  and  time-dependent  sources.  Interestingly,  in  preceding 
sections,  the  approach  consists  of  four-dimensional  quantities,  and  hence  continues  a 
covariant  development. 

5.0  ELECTROMAGNETIC  FIELD  POINT  SPREAD 
FUNCTION  REPRESENTATION 

Inverse-solution  (equation  49)  is  essentially  the  Fourier  transform  of  the  product  of 
two  functions.  By  the  convolution  theorem  (Brace well,  1965)  this  equals  the  convolution 
of  the  Fourier  transforms  of  the  two  functions.  The  convolution  is  then  in  direct-space 
( and  not  K  -  space),  and  hence  is  in  terms  of  source  functions  Sj  ands2  directly.  We 
develop  such  a  direct-space  solution  next. 

We  note  the  following  Fourier  transform  pair  relations 


G(R)  = 


,4k 


+  oo 


1  +.(K.R)(JK 


(-=-5) 

n  K 2 


1 


(~)  =  ( 

7t  K2  4 n 


m 


G  (R)e  nK*R)dR 


for  a  Green's  function  (Jackson,  1962)  or  point  spread  function  given  by 


(50a) 


(50b) 


G  (R)  = 


8(^-0 

c 


(51) 
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The  greek  letter  8  denotes  Dirac  delta  function.  Next  we  form  the  following  partial 
derivatives  of  the  Green's  function  from  equation  50a, 


M 


+  °°  fi  k  \  .  T.  „ 

->  Vn  K2) 

for  u  =  x,y,  z,  orx 


(52a) 


Equation  52a  shows  that  ^-G  ( R )  has  as  its  Fourier  transform  mate  the  quantity 

ou 


m 


^-G(R)e  nK*R)dR. 
ou 


for  u  =  x,y,z,orx 


(52b) 


Note,  quantity  k  in  identity  (equation  52b)  represents  the  general  element  of  matrix 
[K]  given  by  equation  42.  With  repeated  use  of  identity  (equation  52b)  for  each  element 
of  [K]  we  form  the  following  matrix  Fourier  transform  pair 


(_=L)  Ki  K*  =  (JL) 

n  K2  K2Kx  y47i2  J 


NiN2  = 

n2nx 


m 


N,N2  e-HK.R)dR 
No  N, 


_-l  e+i{K.R)dK 

n  K 2  K2  Kx 


(53a) 


(53b) 


where 


n2n 


N,  N, 

..  2  =  IN)  = 


Qx  0  0  -Qx  0  ~QZ  Qy  0 

0  Qx  0  -Qy  Qz  0  -Qx  0 

0  0  Q%-Qz-QyQx  0  0 

Qx  Qy  Qz  Qx  0  0  0  0 

0  -Qz  Qy  0  Q,  0  0  -qx 

Qz  0  -Qx  0  0  Q,  0  -Qy 

-Qy  Qx  0  0  0  0  Qx  -Q. 

o  0  0  0  Q,  Q„  Q,  Qt 
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with  elements 


Qu(R)^<R>- 


The  partial  derivatives  of  the  Green's  function  are  given  by 


tcG(R)  = 


VW3J  c  lc|r|2J  c 


3y  llrl3J  c  clrl2  c 


AG(R)  = 


lw3J  c  Uw2J  c 


j^(R)  =  (J_,s.(y-t). 


(55) 

(56a) 

(56b) 

(56c) 

(56d) 


According  to  plan,  we  now  apply  the  convolution  theorem  to  inversion  formula 
(equation  49)  which  becomes  the  convolution  of  [s^]  T  with  the  Fourier  transform  of 
-  [K]  /  (kK2)  .  By  identity  (equation  53a)  the  latter  is  [N].  Hence  the  inversion  formula 
becomes 


?1<K)1  =  (1)  f  f  f  r+0°k(K-R’)  N2(R-R'j 

[2(r)\  c  j  J  J  J  ~  [n2(r-r')  n^r-r') 


(57) 


Integration  of  equation  57  over  time  yields  the  desired  inversion  formula  in  direct-space. 
The  fields  are  the  results  of  three  added  convolution  effects 


/l(Af)  +  /j  (r,  t)  +  fl(r,t) 

h  0_  f2  (r,  t)  }  f2  (r,  t)  2  +  f2  (r,  t)  ^ 
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where  the  individual  convolutions  are 


A  <r>  o 

f_2  l\  j 


A  <'•  0 


-1 2 


m 

jij 


(r\  f)“ 

*  VI 
N> 

1 _ 

dr’. 


(59a) 


ret 


“I 

dr\ 


(59b) 


ret 


and 


A  (f’ f) 

A  (r’ *) 


-*3 


JJj: 


Z^r-r')  Z2(r-r') 


- 

- 

_^S2  ('•■'•) 

dr. 


(59c) 


Jret 


Quantities  X Y.,  Z i  =1,2  are  defined  by 


Xj  (r) 

X2(r) 

1 

/ 

H 

X2(r) 

X!  (r) 

~  cW3 

H 

/ 

Vj  (r) 

Y2(r)' 

1 

/ 

H 

V'> 

(r) 

-2i  ,i  2 

C  |t’| 

H 

/ 

Zj  (r) 

Z2(r) 

i 

/ 

O 

Z2(r) 

Z  j  (r) 

C2|f| 

0 

/ 

(60a) 


(60b) 


(60c) 
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where 


ro 

0 

0 

*1 

’  0 

z 

-y 

1  = 

j  0 

0 

0 

•ts 

c 

H  = 

-z 

0 

X 

J 

0 

0 

0 

Z  1 

y 

-x 

0 

—x 

-y 

-z 

0 

_  0 

0 

0 

(61a) 


I  = 


10  0  0 
0  10  0 
0  0  10 
0  0  0  1 


O  = 


0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 


(61b) 


The  brackets  []  ret  appearing  in  equation  59  means  that  the  time  t'  is  to  be  evaluated  at 
the  retarded  time 

=  (62) 
c 

Inversion  formula  (equations  58  and  59)  for  the  fields  appears  to  be  new.  We  emphasize 
the  inversion  formula  calculates  the  fields  directly  from  the  sources,  with  no  need  to 
compute  potential  functions  first.  From  the  form  of  equations  59a,  b,  and  c,  quantities 
xfYrzi  play  the  role  of  effective  Green's  functions,  or  point  spread  functions,  for  the 
inversion  process. 

6.0  EXAMPLE:  INFINITE  LINE  OF  ELECTRIC  CHARGE 


Merely  to  exemplify  use  of  the  inversion  approach  (equations  58  through  61),  we 
apply  the  approach  to  the  following  simple  problem.  An  infinite  line  of  electric  charge  is 
uniformly  distributed  along  the  z-axis,  with  a  time-independent  electric  charge  density 
function  given  by 

pV.O  =  P e(r’)  =  X  80t')8(y’).  (63) 

The  constant,  X  ,  represents  the  electric  charge  per  unit  length.  Since  the  charge- 
distribution  is  time-independent,  then  it  is  only  necessary  to  use  the  first  convolution 
integral  in  equation  58  in  calculating  the  electromagnetic  fields. 


fi  (r>  0 

Mr>*) 


III: 


Xj  (r-r’)  X2(r-r’)  Sj(r\f) 

X2(r-r')  X,  (r-r')  s2(r',f) 

1  J  Lz  J  ret 
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The  charge  and  current  density  vector  in  equation  64  for  this  example  is  given  by 


We  can  easily  show  using  equations  65,  60a,  and  61a  that 

+  (x-x'j 

+  (y-y') 

Xl(r-r')  X7(r-r')  s^r'J')  _  tXS  (*’)  5  (y1)  +  2  ) 

X2(r-r')  Xj(r-r')J  |_s2(r’,f’)J  |r-r'J3  0 

ref  Q 

0 

0 


Substitution  of  equation  66  back  into  equation  64  leads  to  the  result 


which  is  equivalent  to  the  well  known  results 

_  ,  v  +  2Xx  _  ,  +2 \y 

E  (r,  f)  =  — 2 2“  =  — 2 2 


U  +y  ] 


Ez(r,t)  =  0 


Hx(r,t)  =  0  H(r,f)=0  Hz(r,t)  =  0. 


(65) 


(66) 


(68a) 

(68b) 
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7.0  ELECTROMAGNETIC  POTENTIAL  FOURIER 
TRANSFORM  REPRESENTATION 


We  showed,  in  sections  4  and  5,  a  solution  for  the  fields  directly  in  terms  of  the 
sources.  Using  the  same  general  approach,  it  is  possible  to  alternatively  express  the 
potential  functions  in  terms  of  the  sources.  Again,  simple  matrix  manipulations  and 
four-dimensional  Fourier  transforms  are  used.  For  brevity,  we  only  give  the  results. 

T* 

Define  [pt  p2]  to  be  the  Fourier  transform  (as  in  equation  38a)  of  the  potentials 
T  T 

[Pj  p2]  •  Then  the  spectra  [Sj  s2]  of  the  sources  connect  with  the  potentials 
through 


Pl 

n 

i* 

~sl 

h 

h 

K2  =  k2+k2+k2-  (c o2/c2) . 
xyz 


(69) 


The  derivation  is  analogous  to  equations  40  through  48.  This  is  a  transfer  theory 
analogous  to  optical-transfer  theory  (Goodman,  1968).  What  might  be  called  an 
"electromagnetic  transfer  function"  is  quantity  l/K  .  As  an  inversion  formula,  equation 
69  would  be  convenient  to  digitally  implement  by  the  use  of  a  four-dimensional  fast 
Fourier  transform  algorithm. 

Solution  (equation  69)  for  the  potentials  corresponds  to  solution  (equation  48)  for  the 
fields.  Both  are  transfer  theorems,  and  both  can  be  ill-posed  as  K  approaches  zero.  A  way 
out  is  to  cast  the  solutions  in  direct-space.  In  the  field  calculations  case,  this  gave  rise  to 
solution  equations  58  and  59.  In  the  potential  calculation  case,  the  corresponding 
solution  is 


Pl  (r,  t) 
P2(M) 


(70) 


This  is  a  well-known  solution  to  the  direct-space  problem. 
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8.0  SUMMARY  AND  CONCLUSIONS 


The  Maxwell  field  equations  have  a  natural  matrix  form  (equation  9)  where  the 
matrix  is  not  square.  With  the  intuitive  goal  of  making  the  matrix  a  square,  and 
continuing  the  skew-Hermitian  trend  of  the  existing  elements  of  the  matrix,  two 
columns  are  added  with  compensating  zeros  in  the  field  vector;  see  equation  10.  Some  of 
the  elements  (see  equations  12a  through  14)  remain  unfixed.  Interestingly,  the  particular 
choices  (equations  13  or  14)  for  the  elements  lead  to  either  the  Lorentz-gauge  or  the 
Coulomb-gauge  representations  of  the  fields,  respectively.  Other  choices  for  the  elements 
will  lead  to  other  gauges. 

The  tactic  of  squaring  matrix  [M]  has  apparently  not  been  done  before,  but  has  many 
benefits,  as  summarized  next.  Use  of  the  square  matrix  operator  [M]  allows  Maxwell's 
equations  to  be  placed  in  a  compact,  covariant  form  (equation  16)  (see  appendix  B  for 
proof  of  covariance).  Once  in  this  form,  use  of  the  simple  matrix  operation  of 
multiplication  upon  equation  16  derives,  in  turn,  the  landmark  effects  of  electromagnetic 
theory:  (a)  the  electromagnetic  wave,  and  charge  continuity  equations,  (equation  20); 
(b)  the  Lorentz-gauge  /  Coulomb-gauge  definitions  (equation  25)  of  the  electromagnetic 
potentials;  (c)  the  wave  equations  (equation  31)  for  the  potentials;  and  (d)  Poyn ting's 
theorem  (equation  34)  on  energy  conservation.  These  matrix-based  derivations 
completely  avoid  the  need  for  vector  calculus  and  attendant  special-purpose  identities 
involving  multiple  curls  and  divergences.  The  necessary  mathematical  baggage  is 
strongly  reduced  by  the  approach.  Note,  the  extension  of  the  matrix  [M]  into  a  square 
form,  described  above,  is  essential  to  the  above  electromagnetic  derivations. 

Alternatively,  taking  the  four-dimensional  Fourier  transforms  of  equation  16  leads  to 
(e)  a  Fourier  representation  (equation  45)  of  Maxwell's  equations;  (f)  their  inversion, 
(equation  48)  for  the  fields  directly  in  terms  of  the  sources  (but  in  Fourier-space);  (g) 
corresponding  inversion  formulae  (equations  58  and  59)  in  direct-space;  and  (h) 
corresponding  results  (see  appendix  A)  to  a  through  c  preceding  for  the  Coulomb  gauge. 

Results  e  through  g  are  new  in  application  to  the  full-fledged  space-time  problem. 
Analogous  results  have  previously  been  found  for  the  special  case  of  time-harmonic 
fields  (Mayes,  1967  and  1985). 
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Other  standard  results  of  electromagnetic  theory  follow  as  well  from  the  matrix 
approach,  but  the  need  for  brevity  rules  out  their  inclusion  here.  For  example,  a  simple 
operation  upon  matrix  equation  16  in  the  source-free  case  directly  leads  to  the  familiar 
requirement  that  the  resulting  electromagnetic  waves  should  be  transverse. 

The  matrix  formulation  involving  the  Lorentz  gauge  is  covariant.  All  basic  vector  and 
matrix  quantities  are  four-dimensional.  However,  the  fact  that  there  are  eight  scalar 
Maxwell  equations  causes  a  requirement  that  all  four-dimensional  quantities  be 
"doubled  up"  into  eight-dimensional  quantities:  Thus  the  two  four-dimensional  sources 
vectors  [sj  ,[s2]  are  stacked  vertically  as  one  eight-dimensional  vector,  see  equation 
16;  the  four-dimensional  matrices  [Mj] ,  [A42]  are  blocked  into  a  single  square  eight-by- 
eight  matrix  [M];  etc.  This  also  emphasizes  that  electromagnetic  theory  is  basically  eight¬ 
dimensional. 

Because  of  the  power  of  matrix  operations,  we  were  able  to  derive  the  key  effects  of 
electromagnetic  theory  without  the  need  for  the  usual  plethora  of  vector  calculus 
identities  that  have  become  the  standard  in  these  derivations:  Stokes'  theorem,  the 
Divergence  theorem.  Green's  theorem,  the  formula  for  the  curl  of  a  curl,  etc.  Instead,  the 
simple  matrix  operations  of  multiplication  and  inversion  were  used.  We  find  these  to  be 
easier  to  understand  and  remember  than  the  preceding  vector  calculus  identities.  Hence, 
we  believe  that  electromagnetic  theory  is  easier  to  learn  (and  is  perhaps  best  taught) 
from  the  matrix  viewpoint. 

Aside  from  pedagogical  aspects,  we  also  found  the  matrix  approach  lends  itself  to 
convenient  linear  analysis  via  Fourier  theory.  Fourier  theory  is,  of  course,  widely  used 
and  understood,  because  of  its  well-known  applications  in  optical  systems  theory 
(Goodman,  1968),  circuit  theory,  etc.  Also,  fast  Fourier  transform  subroutines  widely 
exist,  so  the  approach  is  practical.  Our  application  leads  to  a  Fourier  theory  of 
electromagnetism,  based  on  the  use  of  four-dimensional  Fourier  transforms.  Maxwell's 
equations,  expressed  this  way  in  K-space,  become  a  simple  transfer  theorem  (equation 
45)  with  a  transfer  function  [K].  Very  conveniently,  [K]  has  a  simple  inverse,  and  this 
allows  inversion  (equation  48)  to  take  place.  Moreover,  the  inversion  formulae  equation 
48  or  equations  58  and  59  that  result  express  the  fields  directly  in  terms  of  the  sources. 
There  is  no  need  to  calculate  intermediary  potentials  first. 
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APPENDIX  A 
COULOMB  GAUGE 


As  previously  indicated,  the  conditions  of  the  Coulomb  gauge  are  satisfied  if  the 
quantities  m 44  and  mgg  appearing  in  the  Maxwell  field  equations  10  are  set  equal  to 
zero.  This  is  shown  next.  With  m 44  =  mgg  =  0  ,  matrix  [M]  takes  the  form 
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Then  equation  25  is  explicitly  the  set  of  equations 


E  =  -  V  tpe  -  (-)  ^-Ae  -  Vx  Am 
c  dt 

H  =  -V(pw-  (-)^-Am  +  VxAe 

C  dt 

and 

V*Ae  =  0 


(A2a) 

(A2b) 


(A3a) 


v.aot  =  0 


(A3b) 
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Equations  A2  define  the  electromagnetic  fields  in  terms  of  the  electromagnetic  potentials, 
and  equations  A3  specifically  define  the  Coulomb  gauge  (Jackson,  1962).  This  is  what  we 
set  out  to  show. 


Next,  we  multiply  both  sides  of  equation  16  by  the  complex  conjugate  of  [M]  given 
by  equation  Al.  We  obtain 
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(Compare  with  the  corresponding  Lorentz  gauge  equation  21.)  Again  D2  is  the 
D'Alembertian  operator  and  V2  is  the  Laplacian  operator.  The  single  matrix  equation 
A4  is  equivalent  to  equations  23  through  24  describing  the  electromagnetic  field  wave 
equations  and  the  electric  and  magnetic  charge  continuity  equations.  These  are  the  same 
as  for  the  Ix>rentz  gauge,  as  can  be  easily  verified. 


Finally,  we  substitute  equation  25  back  into  Maxwell's  equations  (equation  16)  with 
[M]  given  by  equation  Al.  We  obtain 
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Pi 
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Pi 

(A6) 


[D]  is  given  in  equation  A5.  This  single  matrix  equation  is  equivalent  to  the  following  set 
of  equations 


W2Ae- =  -  A/  +  (t)V  (Ii<pe) 


c2dt2 


kat 


V2Am_  (  1  jJVi  =  _  (%f>  +  (i)V  (|V”) 

c2  a t2  c  c  dt 


(A7a) 

(A7b) 


and 


V2cp^  =  -(47i)  pe 
V2(pm  =  -(47t)pm. 


(A8a) 

(A8b) 


(Compare  with  the  Lorentz  gauge  equations  32  through  33.)  Equations  (A7)  represent 
the  inhomogeneous  wave  equations  for  the  electric  and  magnetic  vector  potentials,  and 
equations  A8  show  the  scalar  electric  and  magnetic  potentials  each  satisfy  Poisson's 
equation.  These  are  the  known  results  for  the  Coulomb  gauge. 
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APPENDIX  B 

COVARIANCE  OF  THE  EIGHT-BY-EIGHT  MATRIX 
REPRESENTATION  OF  MAXWELL'S  EQUATIONS 


In  this  appendix  we  show  the  eight-by-eight  matrix  representation  (equation  16)  of 
the  Maxwell  field  equations  is  in  covariant  form.  Consider  two  inertial  right-handed 
Cartesian  coordinate  reference  frames  S  and  S'  as  shown  in  the  figure  below.  These  two 
frames  share  a  common  zz'  axis.  The  frame  S'  is  receding  from  the  frame  S  with  a 
constant  speed  v  in  the  positive  z-direction. 


y  y' 


Figure  Bl.  Two  inertial  reference  frames  S  and  S'  in  relative  motion. 

An  observer  in  S  reports  space-time  coordinates  (x,y,z,t)  for  an  event  and  an  observer  in 
S'  reports  (x  ,y' ,z  ,?)  for  the  same  event.  For  convenience  we  suppose  that  each  observer 
chooses  t  =  t'  =  0  to  represent  the  instant  that  their  origins  coincide.  We  know  that  the 
space-time  coordinates  in  S  and  S'  are  related  through  the  Lorentz  transformation 
(Jackson,  1962). 


4 

x  jx  =  ^  ^pv*v'  ^  =  L2, 3, 4  (Bl) 

v  =  l 

with  relativistic  four-vectors 


(xvx2,xvx4)  =  (x,y,  z,x),  ( x\,x'2,x\,x\ )  =  (x\  y\  z\  x’) . 


(B2) 
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Coefficients  L 


form  the  four-by-four  Lorentz  transformation  matrix 
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The  speed  parameter,  (3  ,  and  the  Lorentz  factor,  y,  are  given  by 

o  v  1 

P  =  7-  y  *  -J=5-  (B4) 

Jl-P2 

The  transformation  properties  of  the  electromagnetic  field  (Jackson,  1962)  can  be  found 
from 


4  4 

^pv  ~  M^>v  =  1>2, 3,4  (B5) 

A,  =  l  o  =  l 

where  are  elements  of  the  anti-symmetric  field-strength  tensor 
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Now,  by  equation  B6  expressed  in  the  primed  system,  element  F’4l  is  iE  ,  f"42  is 

'l 

iEy  ,  etc.  Note,  that  elements  /£'  ,  etc.,  define  vector  [f\  /^l  .  Hence, 
transformation  equation  B5  may  be  used  to  link  [f\  f2] T  to  a  linear  superposition  of 
unprimed  elements  F^a  which,  again,  by  equation  B6  are  elements  of  \fx  f2) T  ■  These 
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are  linear  equations  of  the  matrix  form 
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It  is  of  interest  to  note  that  the  length  of  the  eight-by-one  electromagnetic  field  vector 
is  invariant  under  the  Lorentz  transformation  (equation  B7).  That  is. 
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For  this  reason,  this  vector  could  be  referred  to  as  an  electromagnetic  relativistic  eight- 
vector.  We  are  not  aware  that  this  concept  has  been  used  before. 

Since  and  are  relativistic  four-vectors,  then  the  eight-by-one  source  vector 
appearing  in  equation  16  transforms  between  the  frames  S  and  S'  according  to 
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where  "element"  O  is  a  four-by-four  matrix  of  zeros  defined  by  equation  21. 

We  are  now  in  a  position  to  show  that  equation  16  is  in  covariant  form.  Starting  with 
the  matrix  representation  of  the  Maxwell  field  equations  in  inertial  frame  S,  namely 
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equation  16,  we  perform  the  following  matrix  multiplication  operation 
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With  the  use  of  the  left-hand  equation  in  equation  BIO,  we  can  rewrite  equation  Bll  in 
the  form 
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Next  we  replace  the  electromagnetic  relativistic  eight-vector  in  equation  B12  by  the  right- 
hand  equation  in  B7.  This  gives 
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With  the  use  of  the  following  partial  differential  operator  transformation  (Jackson,  1962) 
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it  can  be  shown  that 
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With  the  use  of  equation  B15,  equation  B13  reduces  to 
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(B15) 


(B16) 


which  we  recognize  as  the  matrix  representation  (equation  16)  of  the  Maxwell  field 
equations  in  frame  S'.  This  establishes  the  matrix  representation  (equation  16)  of  the 
Maxwell  field  equations  is  in  covariant  form. 
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APPENDIX  C 

INVERSION  OF  MATRIX  [K] 


Next,  consider  the  left-hand  side  of  equation  Cl.  Matrix  [M]  follows  the  ordinary 
matrix  rules  of  multiplication,  so  that  it  obeys  the  associative  law.  Then  the  left-hand  side 
of  equation  Cl  obeys 
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Now  the  product  [M]*[M]  obeys,  by  equations  20  and  21 
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Hence,  the  left-hand  side  obeys. 
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By  straightforward  differentiation  and  the  use  of  notation  in  equation  61b 
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where  K 2  is  the  simple  scalar  given  by  equation  47.  Then  combining  equations  C8  and 
C9,  the  left-hand  side  is 
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Then  since  the  left-hand  side  and  right-hand  sides  of  equation  Cl  are  equal,  by 
equations  C5  and  CIO 
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By  definition  of  the  inverse,  this  implies  that 


,-2 

*1  *2 
*2  *1 


K2 


*1  *2 
*2*1 


(Cl  3) 


C-3 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  Information  Is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information.  Including 
suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services.  Directorate  for  Information  Operations  and  Reports.  1215Jefferson  Davis  Highway.  Suite  1204,  Arlington.  VA  22202-4302. 
and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188).  Washington.  ‘:»C  20503  _ _____ _ _ 


2  REPORT  DATE  3  REPORT  TYPE  AND  DATES  COVERED 


1  AGENCY  USE  ONLY  (Leave  blank) 


March  1992 


4  TITLE  AND  SUBTITLE 

A  NEW  MATRIX  FORMULATION  OF  CLASSICAL  ELECTRODYNAMICS 
Part  1.  Vacuum 


6  AUTHOR(S) 

R.  P.  Bocker  and  B.  R.  Frieden 


6  FUNDING  NUMBERS 


In-house 


7  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Command,  Control  and  Ocean  Surveillance 
Center  (NCCOSC),  Research,  Development,  Test 
and  Evaluation  Division  (NRaD) 

San  Diego,  CA  92152-5000 


9  SPONSORING/MONITORING  AGENCY  NAME<S)  AND  ADDRESS(ES) 


University  of  Arizona 
Optical  Sciences  Center 
Tucson,  Arizona  85721 


8  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

NRaD  TD  2259 


10  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 


12a.  DISTRIBUTION/AV A! LABILITY  STATEMENT 


Approved  for  public  release;  distribution  is  unlimited. 


13  ABSTRACT  (Maximum  200  words) 

Presented  in  this  paper  is  a  new  matrix  representation  of  classical  electromagnetic  theory.  The  basis  of  this  representa¬ 
tion  is  a  space-time,  eight-by-eight  differential  matrix  operator.  This  matrix  operator  is  initially  formulated  from  the  differ¬ 
ential  form  of  the  Maxwell  field  equations  in  vacuum.  The  resulting  matrix  formulation  of  Maxwell’s  equations  allows  sim¬ 
ple  and  direct  derivation  of  the  electromagnetic  wave  and  charge  continuity  equations,  the  Lorentz  conditions  and  defini¬ 
tion  of  the  electromagnetic  potentials,  the  Lorentz  and  Coulomb  gauges,  the  electromagnetic  potential  wave  equations,  and 
Poynting’s  conservation  of  energy  theorem.  A  four-dimensional  Fourier  transform  of  the  matrix  equations  casts  them  into 
an  eight-dimensional  transfer  theorem.  The  transfer  function  has  an  inverse,  and  this  allows  the  equations  to  be  inverted. 
This  inversion  expresses  the  fields  directly  in  terms  of  the  charge  and  current  source  distributions,  i.  e.,  without  the  need 
for  calculating  intermediary  potentials.  This  inversion  formula  is  new,  for  the  general  scenario  of  space-  and 
time-dependent  sources.  A  simple  pedagogical  example  is  included  illustrating  use  of  the  formulation. 


17  SECURITY  CLASSIFICATION 
OF  REPORT 


18  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 


19  SECURITY  CLASSIFICATION 
OF  ABSTRACT 


15  NUMBER  OF  PAGES 


16  PRICE  CODE 


20  LIMITATION  OF  ABSTRACT 


UNCLASSIFIED 


UNCLASSIFIED 


UNCLASSIFIED 


SAME  AS  REPORT 


NSN  754001  280-5500 


Standard  form  298 


UNCLASSIFIED 


NSN  75*0-01-280-5600 


Standard  toon  ZM 

UNCLASSIHED 


INITIAL  DISTRIBUTION 


Code 

0012 

Patent  Counsel 

(1) 

Code 

0141 

A.  Gordon 

(1) 

Code 

414 

T.  Phillips 

(1) 

Code 

635 

H.  E.  VanBrundt 

(1) 

Code 

653 

K.  Swecker 

(1) 

Code 

70 

E.  Shutters 

(1) 

Code 

743 

H.  J,  Whitehouse 

(1) 

Code 

75 

J,  E.  Griffin 

(1) 

Code 

7502 

B.  R.  Hunt 

(1) 

Code 

754 

D.  K.  Forbes 

(1) 

Code 

76 

J.  R.  Wangler 

(1) 

Code 

7601 

J.  R.  Zeidler 

(1) 

Code 

7601 

K.  Bromley 

(1) 

Code 

7602 

S.  Hoff 

(1) 

Code 

761 

H.  M.  Berry 

(1) 

Code 

761 

R.  P.  Bocker 

(60) 

Code 

761 

G.  W.  Byram 

(1) 

Code 

761 

S.  I.  Chou 

(1) 

Code 

761 

M.  Lasher 

(1) 

Code 

761 

J.  P.  Loughlin 

(1) 

Code 

765 

E.  Munn 

(1) 

Code 

80 

K.  D.  Regan 

(1) 

Code 

805 

J.  W.  Rockway 

(1) 

Code 

808 

E.  Schimitschek 

(1) 

Code 

808 

G.  W.  Beaghler 

(1) 

Code 

82 

R.  J.  Ko chan ski 

(1) 

Code 

824 

J.  C.  Logan 

(1) 

Code 

832 

P.  M.  Hansen 

(1) 

Code 

84 

M.  S.  Kvigne 

(1) 

Code 

844 

B.  James 

(1) 

Defense 

Technical  Information  Center 

Alexandria,  VA  22304-6145  (4) 


NCCOSC  Washington  Liaison  Office 
Washington,  DC  20363-5100 

Center  for  Naval  Analyses 
Alexandria,  VA  22302-0268 

Navy  Acquisition,  Research  &  Development 
Information  Center  (NARDIC) 

Alexandri,  VA  22333 

Navy  Acquisition,  Research  &  Development 
Information  Center  (NARDIC) 

Pasadena,  CA  91106-3955 


